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Abstract 


The distribution and orientation of energy inside jets is predicted to be an experimental 
handle on colour connections between the hard-scatter quarks and gluons initiating the jets. 
This Letter presents a measurement of the distribution of one such variable, the jet pull 
angle. The pull angle is measured for jets produced in tt events with one W boson decaying 
leptonically and the other decaying to jets using 20.3 fb” 1 of data recorded with the ATLAS 
detector at a centre-of-mass energy of yfs = 8 TeV at the LHC. The jet pull angle distribution 
is corrected for detector resolution and acceptance effects and is compared to various models. 
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1 Introduction 


Due to the confining nature of the strong force, directly measuring the quantum chromodynamic (QCD) 
interactions between quarks and gluons is not possible. The strength and direction of the strong force 
depends on the colour charge of the particles involved. To a good approximation, the radiation pattern in 
QCD can be described through a colour-connection picture, which consists of colour strings connecting 
quarks and gluons of one colour to quarks and gluons of the corresponding anti-colour. An important 
question is whether there is evidence of these colour connections (colour flow) in the observable objects: 
colour-neutral hadrons and the jets they form. The study of energy distributions inside and between 
jets in various topologies has a long history, dating back to the discovery of gluons in three-jet events 
at PETRA [1-4]. Colour connections are still a poorly constrained QCD effect, which motivates the 
dedicated study presented in this Letter. If well understood, experiments can exploit colour flow to aid 
Standard Model (SM) measurements and searches for physics beyond the SM (BSM). As a test that 
the colour flow can be extracted from the observable final state, the data are compared to models with 
simulated W bosons that are colour-charged or colour-neutral. 


One observable predicted to contain information about the colour representation of a dijet resonance like 
the W, Z, or Higgs boson, is the jet pull vector [5]. The pull vector for a given jet J with transverse 
momentum, is defined as 
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The sum in Eq. (1) runs over jet constituents with transverse momentum p' T and location r\ = (Ay/,, A<j>j), 
defined as the vector difference between the constituent and the jet axis (yj,(f>j) in rapidity ( y ) - azimuthal 
angle (c/>) space 1 . Given the pull vector for jet J\, the angle formed between this pull vector and the vector 
connecting J\ and another jet T, (yj 2 - yj ] , cbj, - (t>j ] ), is expected to be sensitive to the underlying colour 
connections between the jets. This is shown graphically in Fig. 1, and the angle is called the pull angle, 
denoted 6p(J \, Jj)- The pull angle is symmetric around zero when it takes values between -n and n and 
so henceforth Op(J\, Jfl) refers to the magnitude of the angle in (Ay, Ac/)) space with 0 < Op < n. If the 
pull vector is computed using jets originating from colour-connected quarks, 6p ~ 0 since the radiation 
is predicted to fall mostly between the two jets. If the vector is computed using two jets which do not 
originate from colour-connected quarks, there is no reason to expect dp to be small. Thus dp should be 
useful for determining colour connections. 


One of the challenges in studying colour flow is the selection of a final state with a known colour compos¬ 
ition. Colour-singlet W bosons from tt events provide an excellent testing ground because these bosons 
have a known initial (colourless) state and such events can be selected with high purity. The first study 
of colour flow using W bosons from top quark decays was performed by the D0 Collaboration [6]. In 
the D0 analysis, calorimeter cells clustered within jets were used as the constituents in Eq. (1) and the 
resulting distribution was compared to W singlet (nominal) and W octet templates. The impact of the 
colour flow on the observable energy distributions is very subtle; the D0 result was statistically limited 
and was not able to significantly distinguish singlet and octet radiation patterns. 


1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the 
detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the i/-axis 

points upward. Cylindrical coordinates (r, <f>) are used in the transverse plane, (j) being the azimuthal angle around the beam 
pipe. The pseudorapidity is defined in terms of the polar angle 6 as 77 = -lntan(0/2). Separation between objects in 
space is defined as A R = -y/f A //) 2 + (A c/>) 2 . The rapidity of a four-vector y = 0.5 log(|^), where E is the energy and p : the 
component of the momentum parallel to the beam axis. 
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Fig. 1: Diagram illustrating the construction of the jet pull angle for jet J\ with respect to Jo- 


The analysis presented in this Letter is a measurement of the jet pull angle in Vs = 8 TeV pp collisions 
at the Large Hadron Collider (LHC) with the ATLAS detector. Instead of comparing the reconstructed 
jet pull angle in data directly to simulated templates, as was done by the D0Collaboration, the jet pull 
angle distribution is first corrected for detector resolution and acceptance effects. This allows for direct 
comparison with particle-level predictions for models of physics beyond the SM as well as simulations 
of non-perturbative physics effects with various tunable parameter values. 


2 Object and event selection 

The ATLAS detector independently measures the inclusive and charged-particle energy distributions in 
jets. This allows the jet pull angle to be constructed using only the charged constituents of jets, or both 
the charged and neutral constituents. In this analysis, both options are used in order to provide inde¬ 
pendent measurements with different experimental systematic uncertainties. Charged-particle momenta 
are measured in a series of tracking detectors (collectively called the inner detector), covering a range of 
|//| < 2.5 and immersed in a 2 T magnetic field. Electromagnetic and hadronic calorimeters surround the 
inner detector, with forward calorimeters allowing electromagnetic and hadronic energy measurements 
up to I 77 I =4.5. A detailed description of the ATLAS detector can be found elsewhere [7]. 

The anti-k r algorithm [ 8 ] with radius parameter R = 0.4 is used to reconstruct jets from clusters of calor¬ 
imeter cells [9] with deposited energy. The clusters are formed using the local cluster weighting (LCW) 
algorithm [ 10 ] and calibrated to account for the detector response as well as to mitigate the contributions 
from additional pp collisions per bunch crossing (pileup) [11]. The all-particles pull angle is built from 
the clusters of calorimeter cells assigned to a given jet. In order to improve the rapidity resolution, jets and 
clusters are corrected to point toward the reconstructed primary vertex 2 and the corrected jet four-vector 
is used as the axis in Eq. (1). 


2 The primary vertex is defined as the vertex with the highest 2 p\ of associated tracks. 
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Process 

Number of events 

tt 

95400 ± 14000 

IT/-channel single top 

2730 ± 600 

s- and (-channel single top 

150 ± 10 

IT+jets 

3710 ± 120 

Z+jets 

560 ± 270 

Dibosons 

190 ± 40 

Multijets 

2500 ±910 

Total SM 

105000 ± 14000 

Data 

102987 


Table 1: Estimated composition of the selected event sample. The uncertainties are the sum in quadrature of the 
statistical uncertainties and either the uncertainties of the normalisation method (for the data driven W+jet and 
multi-jet estimates) or the uncertainties of the cross-section estimates. 


The charged-particles pull angle is built from tracks that are associated with a given jet [12, 13]. Tracks 
are reconstructed from hits in the inner detector and are required to have |//| < 2.5, p-y > 0.5 GeV, and 
satisfy various quality criteria, such as association with the primary vertex, in order to suppress tracks 
from random hits and pileup tracks. The charged-particles pull angle is constructed using the four- 
momentum sum of all the associated tracks (treated as massless) to provide the axis in Eq. (1). The 
all-particles and charged-particles pull angles are largely independent as they rely on information from 
different detector sub-systems. 

In order to isolate a pure sample of hadronically decaying W bosons, this analysis targets a it —» bbW(—> 
qq')W(—> tv) final state. 

Events arc selected by triggers requiring a single isolated electron or muon, and the offline ‘tight’ elec¬ 
tron [14] or ‘combined’ muon [15] must have p-y > 25 GeV and |//| < 2.5. Basic quality criteria are 
imposed, including the existence of at least one primary vertex associated with at least five tracks with 
Pt > 0.4 GeV. Furthermore, the magnitude £™ ss of the missing transverse momentum [16] is required 
to be greater than 20 GeV and Zi™ ss + m T > 60 GeV, where m-y is the transverse mass of the selected 

lepton and the i?™ ss , m T = ^2p e T Ej lss (l - cos(0 f - cp v )). Events must have > 4 jets with px > 25 GeV. 
At least two of these jets must be tagged as 7>-jets using the multivariate discriminant, MV1 [17], which 
uses impact parameter and secondary vertex information. The chosen MV 1 working point corresponds 
to an average ^-tagging efficiency of 70% for A jcts in simulated it events. At least two jets must not be 
^-tagged; of these, the two leading-px jets with \q\ <2.1 are labelled as the jets from the hadronically 
decaying W boson,./, with p^ > p^ . The ^-tagged jets and the jets selected for the pull angle calculation 
are required to have |//| < 2.1 so that all constituents are within coverage of the inner detector used for 
tracking. This procedure selects a sample that is expected to contain approximately 90% tt events. In 
45% of simulated tt events both jets selected for the pull angle calculation contain energy from the had¬ 
ronically decaying W boson. Table 1 shows the predicted composition compared to the data yield. More 
details about the various contributions are given below. 

There are two pull angles calculated per event, one calculated from the all-particles pull vector and one 
from the charged-particles pull vector of the highest pj jet assigned to the hadronic W boson decay. The 
all-particles and charged-particles pull angles are the angles that the corresponding pull vectors make 
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Fig. 2: The detector-level (a) all-particles and (b) charged-particles pull angle, dp, in data and in simulation. 
The uncertainty band includes only the experimental uncertainties on the inputs to the event selection and the jet 
pull calculation. A large part of the uncertainty displayed here affects the overall normalisation and is correlated 
between the individual bins. This component of the uncertainty is cancelled in the unfolded measurement of the 
unit-normalised pull angle distribution. 


with the direction from J\ to J 2 . Figs. 2(a) and 2(b) show comparisons between data and simulation 
for the all-particles and charged-particles pull angles, calculated at detector level, i.e. from reconstructed 
objects. Both distributions are broadly flat with an enhancement at small angles, consistent with the SM 
prediction. 


3 Event simulation 

Monte Carlo (MC) samples are produced in order to determine how the detector response affects the pull 
angle and to estimate some of the non-ff contributions in the data. The details of the samples used are 
shown in Table 2. 

Aside from the VT+jcts background, all MC samples are normalised to their theoretical cross-sections, 
calculated to at least next-to-leading order (NLO) precision in QCD [33-38]. For the purpose of com¬ 
parison between data and the SM prediction before unfolding, tt events are normalised to a cross-section 
of 253 ± 15 pb, calculated at next-to-next-to-leading order (NNLO) in QCD including next-to-next-to- 
leading logarithmic (NNLL) soft gluon terms [39], assuming a top-quark mass of 172.5 GeV. 

Generated events are processed with a full ATLAS detector and trigger simulation [40] based on Geant4 [41] 
and reconstructed using the same software as the experimental data. The effects of pileup are modelled by 
adding to the generated hard-scatter events multiple minimum-bias events simulated with Pythia 8.160 
[42], the A2 set of tuned MC parameters (tune) [43] and the MSTW2008LO Parton Distribution Func¬ 
tion (PDF) set [44]. The distribution of the number of interactions is then weighted to reflect the pileup 
distribution in the data. 
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Process 

Generator 

Type 

Version 

PDF 

Tune 

tt 

Powheg [18-20] 

NLO ME 


CT10 [21, 22] 


+PYTHIA [23] 

+ PS 

6.426.2 

CTEQ6L1 [24] 

Perugia2011c [25] 

Single top 

Powheg 

NLO ME 


CT10(4f) 


+Pythia 

+ PS 

6.426.2 

CTEQ6L1 

Perugia 201 lc 

WW, WZ, ZZ 

Sherpa [26] 

LO multi-leg ME + PS 

1.4.1 

CT10 

Default 

W/Z+jets 

Alpgen [27] 

LO multi-leg ME 

2.1.4 

CTEQ6L1 

- 

+PYTHIA 

+ PS 

6.426.2 

CTEQ6L1 

Perugia 201 lc 


Powheg 

NLO ME 


CT10 


tt f 

+Herwig [28] 

+ PS 

6.520.2 

CT10 

AUET2 [29] 


+Jimmy [30] 

(MPI) 

4.31 

- 

- 


MC@NLO[31, 32] 

NLO ME 

4.06 

CT10 

- 

tt f 

+Herwig 

+ PS 

6.520.2 

CT10 

AUET2 


+Jimmy 

(MPI) 

4.31 

“ 

“ 


Table 2: Monte Carlo samples used in this analysis. The abbreviations ME, PS, PDF, MPI, LO and NLO respectively 
stand for matrix element, parton shower, parton distribution function, multiple parton interactions, leading order and 
next-to-leading order in QCD. Tune refers to the used set of tunable MC parameters. Those samples marked with 
a f are used as alternative tt samples to evaluate uncertainties due to the modelling of tt events. The nominal tt 
generator in the first line is used to estimate the yield in Table 1. 


To test the sensitivity of the jet pull angle to the singlet nature of the W boson, a sample was generated 
with a colour-octet W boson. The octet W boson is simulated using the same setup as the nominal tt setup 
described in Table 2. Using the partons produced with Powheg (in the LHE [45] format), the colour flow 
is inverted such that one of the W decay daughters shares a colour line with the ft-quark and the other 
shares a line with the top quark, as demonstrated schematically in Fig. 3. This sample is referred to as 
colour flipped in the rest of this Letter. 


4 Unfolding 


In order to make direct comparisons with various theoretical models, the data are unfolded to particle- 
level objects. Particle-level jets are clustered from simulated particles with a mean lifetime t > 30 ps, 
before taking into account interactions with the detector. The particle-level inputs to the all-particles pull 
angle arc all of the charged and neutral particles clustered within particle-level jets. Only the charged 
particles clustered within the particle-level jets arc used for the charged-particles pull angle. Additional 
information about the particle-level object definitions can be found in Ref. [46]. 

The particle-level event selection is analogous to the detector-level selection described in Sec. 2 with 
detector-level objects replaced with particle-level objects. Exactly one electron or muon and at least four 
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Fig. 3: Diagram illustrating the colour connections (thick coloured lines) for the nominal sample with a colour- 
singlet W (left) and the flipped sample with a colour-octet W (right) samples in cases where the colour is distorted. 


jets are required, each with p-y > 25 GeV and |//| < 2.5, with events discarded if the electron or muon is 
within A R = 0.4 of a jet. The particle-level missing transverse momentum, E,™ ss , calculated using the 
sum of neutrino four-momenta, is required to be £'™ 1SS > 20 GeV and the sum of is™ ss + mj > 60 GeV. At 
least two of the selected jets are required to be identified as Z?-jets using the same definition as that found 
in Ref. [46]. As with the detector-level calculation of the pull angle, the two leading-px particle-level 
non b -jets with |//| < 2.1 are labelled as the jets from the hadronically decaying W boson. 

The first step of the unfolding procedure is to subtract from the data an estimate of the non-tt backgrounds 
bin-by-bin in the pull angle distribution. The VV'+jcts and multijet backgrounds are estimated from the 
data using the charge asymmetry and matrix methods, respectively [47]. The other backgrounds are estim¬ 
ated from simulation. Single-top Wt events have hadronically decaying W bosons and are thus sensitive 
to colour flow; however, even large changes in the colour flow through the single-top contribution are 
subdominant compared to other uncertainties, and are thus ignored. 

After subtraction, the data are unfolded using an iterative Bayesian (IB) technique [48] as implemented 
in the RooUnfold framework [49]. The IB method iteratively applies Bayes’ theorem using the response 
matrix to connect the prior to posterior at each step. The response matrix is constructed from the nominal 
tt simulation (also used for the prior) and describes the bin migrations between the particle-level and 
detector-level pull angle distributions. The binning of the response matrix and the number of iterations 
in the IB method are chosen to minimise the overall uncertainty, described in Section 5. In addition 
to correcting for bin migrations, the unfolding procedure also corrects for events that pass either the 
detector-level or particle-level selection, but not both. The corrections are estimated using simulated tt 
events. Approximately 70% of events that pass the detector-level selection also pass the particle-level 
selection, while approximately 20% of events that pass the particle-level selection also pass the detector- 
level selection. These corrections are found to be largely independent of the pull angle. 

The all-particles pull angle has a stronger dependence on the colour flow, but has a worse resolution than 
the charged-particles pull angle 3 and is more sensitive to the particle-level spectrum used as the IB prior. 
Three bins and 15 iterations are used for the all-particles pull angle and four bins and three iterations 
are used for the charged-particles pull angle. More iterations are required for the all-particles pull angle 
to reduce the sensitivity to the IB prior, at the cost of a larger statistical uncertainty. The size of each 
bin is comparable to the pull angle resolution. The uncertainty on the unfolding procedure is determined 


3 The RMS of the all-particles (charged-particles) pull angle is found to be 0.35;r (0.287T) radians. 
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in a data-driven way by reweighting the particle-level distributions so that the corresponding detector- 
level distributions are in better agreement with the data [50]. The method uncertainty is then estimated 
using the nominal response matrix by comparing the unfolded, reweighted simulation with the reweighted 
particle-level distribution. 


5 Systematic uncertainties 


There are various systematic uncertainties which can impact the measured jet pull angle. The sources 
of uncertainty can be classified into two categories: experimental uncertainties and theoretical modelling 
uncertainties. In the first category, some uncertainties impact the pull angle directly while the others 
impact only the acceptance. 

The uncertainty in the energy scale and angular distribution of clusters of calorimeter cells are primary 
sources of uncertainty that directly impact the pull angle. To estimate the uncertainty in the angular resol¬ 
ution of clusters, isolated tracks are matched to clusters in events enriched in Z(—> pp) + jets events. The 
distribution of AR( track, cluster) is investigated for various bins of track p-\ and //, and in extrapolations 
to various layers of the calorimeter. The observed resolution uncertainty is < 1 mrad, but to account for 
the possible effects of multi-particle clusters, 5 mrad is used, as in Ref. [51]. Estimations of the cluster 
energy scale uncertainty are based on measurements of the ratio of the cluster energy to the pt of matched 
tracks (E/p) [52]. The ratio E/p in 8 TeV data is compared with simulation and the differences are used to 
estimate the uncertainty, which is parameterised as a( I +/i/p), where a and fi depend on //. Typical values 
are a = 0.05 and ft = 0.5 GeV . To account for an uncertainty in the energy lost due to inactive material 
in the detector and noise thresholds, low-energy clusters arc randomly removed with energy-dependent 
probabilities, as in Ref. [53], 

For the charged-particles pull angle, there are uncertainties associated with track reconstruction effi¬ 
ciency, which are estimated by randomly removing tracks with an //-dependent probability [54]. Fur¬ 
thermore, tracks are randomly removed from jets with a px-dependent factor to estimate the effect of 
uncertainties in the efficiency of track reconstruction inside jets. For jets with pj < 500 GeV, the uncer¬ 
tainty on the track reconstruction efficiency is less than 1%. 

As the pull vector definition uses the calorimeter jet pt, both the all-particles and charged-particles 
pull angle are affected by the uncertainty in the jet energy scale [55, 56] and resolution [57]. However, 
changes in the jet energy scale and resolution do not impact the pull angle, but do impact the results via 
the acceptance due to pj thresholds. Similarly, uncertainties in the lepton energy scale, trigger efficiency, 
Zi™ ss resolution and //-tagging efficiencies [17, 58-60] indirectly affect the results through changes in 
acceptance. Other sources of uncertainty on the acceptance, which impact the measurement through 
the background subtraction, include those related to the luminosity [61], the multijet estimation, and the 
normalisation and heavy flavour content of the VT+jcts background [47]. 

The modelling of the tt sample is a primary source of theoretical uncertainty. For instance, there are 
different treatments of the matrix element calculation (Powheg+Herwig versus MC@NFO+Herwig) and 
the parton shower model (Powheg+Pythia versus Powheg+Herwig). The change in the result when us¬ 
ing a response matrix constructed from the flipped tt simulation is considered as a source of uncertainty 
on the colour model. Other sources of modelling uncertainties are evaluated, including initial/final-state 
radiation by varying the radiation simulated with AcerMC 3.8 [62] constrained by Ref. [63], multiple 
simultaneous parton collisions using the Perugia 2012 mpiHi tune [25], colour reconnections with the 


Source 

All particles [%] 

Charged particles [%] 

Unfolding procedure 

0.5 

0.6 

Clusters 

0.5 

N/A 

Tracks 

N/A 

0.2 

JES 

0.4 

0.2 

JER 

0.3 

0.1 

Matrix element 

1.5 

0.9 

Shower model 

1.7 

0.6 

Colour model 

1.3 

1.0 

ISR/FSR 

1.2 

0.2 

MPI 

0.1 

0.6 

Other 

0.4 

0.2 

Total systematic uncertainty 

3.0 

1.8 

Statistical uncertainty 

1.1 

0.7 


Table 3: Summary of systematic uncertainties in the first bin of the all-particles and charged-particles pull angle 
distributions. The jet energy scale, jet energy resolution, multiple simultaneous parton collision modeling, and 
initial/final-state radiation modeling uncertainties are abbreviated as JES, JER, MPI, and ISR/FSR, respectively. 


Perugia 2012 lowCR tune, overlap between single top and tt by comparing the diagram removal scheme 
with the diagram subtraction scheme [64], and PDF uncertainties [65]. Colour flow describes the colour 
representation of the hard-scatter process while colour reconnection is a phenomenological model for 
implementing a finite number of colours in the parton shower. The observable consequences of both 
processes are similar, but colour reconnection has been constrained by charged-particle distributions in 
experimental data [25]. Varying the top quark mass by ±1 GeV has a negligible impact on this measure¬ 
ment. 

The systematic uncertainties are estimated by unfolding the data with varied response matrices or by 
subtracting varied background predictions from the data. The sources of uncertainty described above are 
summarised in Table 3 for the first bin of the pull angle for both the all-particles and charged-particles 
pull angle. In general uncertainties are found to be larger for the all-particles pull angle than the charged- 
particles pull angle due to the remaining sensitivity to the IB prior, as discussed in Section 4. All bins 
carry information about the underlying colour flow, but the first bin has the largest expected difference 
between the flipped and nominal models. 


6 Results 

The unfolded data are shown in Fig. 4 (all-particles) and Fig. 5 (charged-particles) and compared to SM 
tt predictions simulated with Powheg+Pythia and Powheg+Herwig, and a flipped tt simulation generated 
using Powheg+Pythia. 

The data favour the SM predictions over the prediction of the flipped model. Of the two SM predictions, 
Powheg+Pythia offers a slightly better description of the data. The level of agreement between the data 
and the Powheg+Pythia models is quantified using Ax 2 as a test statistic. A /j-valuc is computed taking 
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into account systematic and statistical uncertainties and their correlations by generating pseudo-datasets 
from the full covariance matrix and then normalising the pull angle distributions. This p-value is then 
converted into a one-sided Gaussian equivalent Z-score. The unfolded data are found to agree with the 
nominal SM colour flow (with Powheg+Pythia) at the 0.8cr (0.9o-) level and differ from the flipped model 
at 2.9<t (3.7ct) for the all-particles (charged-particles) pull angle. Although the statistical uncertainty is 
largely uncorrelated between the all-particles and charged-particles pull angle measurements, the sys¬ 
tematic uncertainty is not. Since the charged-particles pull angle measurement is more sensitive than 
the all-particles pull angle measurement, a full combination of the two measurements does not provide 
increased sensitivity. 


7 Summary 


The analysis presented in this Letter describes a measurement of the orientation of radiation from jets 
identified as originating from a W boson in tt events. The measurement uses 20.3 ftr 1 of yfs = 8 TeV 
pp collision data recorded by the ATLAS detector at the LHC. To quantify the distribution of energy 
inside one jet relative to another, the distribution of the jet pull angle is extracted from the data using 
information from both the ATLAS calorimeter and tracking detectors. The jet pull angle is found to 
correctly characterise the W boson as a colour singlet, with data disfavouring an alternative colour-octet 
model at greater than 3<x. This illustrates the potential to use the jet pull angle in future SM measurements 
and BSM searches. The jet pull angle measurement is presented as a normalised fiducial tt differential 
cross-section, allowing the results to be used to constrain implementations of colour connection. 
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Fig. 4: Normalised fiducial tt differential cross-section for the jet pull angle distribution constructed using all 
particles. The data are compared to two models with a colour-singlet W boson (SM) but with different MC gener¬ 
ators, and to the flipped model with a colour-octet W boson. 
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Fig. 5: Normalised fiducial tt differential cross-section for the jet pull angle distribution constructed using charged 
particles. The data are compared to two models with a colour-singlet W boson (SM) but with different MC gener¬ 
ators, and to the flipped model with a colour-octet W boson. 
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